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SUMMARY 

v 

Helicopter  flight  data,  especially  measurements  of  linear 
accelerations  and  angular  velocities,  are  typically  corrupted  by  sinusoidal 
deterministic  disturbances,  that  are  associated  with  the  rotor  frequency  and 
its  harmonics.  In  an  effort  to  eliminate  these  disturbances  without  distorting 
the  underlying  trends,  a  new  digital  notch  filter,  developed  under  a  research 
agreement  between  Aeronautical  Research  Laboratories  (ARLV4  and  the 
University  of  Newcastle,  has  been  modified  and  implemented  on  the  ARL 
ELXSI  6400  computer.  In  this  memo,  the  filter  is  described  and  the  frequency 
and  transient  response  characteristics  are  summarised.  Practical 
considerations  arising  out  of  application  of  single  and  cascaded  filters  to  Sea 
King  flight  data  are  discussed,  and  the  performance  of  the  new  filter  is 
compared  with  that  of  comparable  Butterworth  filters. 
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1 .  INTRODUCTION 

In  1979,  flight  trials  of  the  Sea  King  Mk  50  helicopter  were 
performed  in  order  to  validate  a  mathematical  model  of  the 
helicopter,  developed  in  the  Aircraft  Behaviour  Studies  -  Rotary 
Wing  (ABS-RW)  Group,  at  Aeronautical  Research  Laboratories  (ARL) 
In  an  attempt  to  remove  unwanted  noise  from  the  data  records,  the 
recorded  signals  were  pre-processed  using  analogue  Butterworth 
filters.  However,  considerable  noise  still  existed  on  the  data 
obtained.  To  remove  this  noise,  one  of  two  digital  Butterworth 
filters  was  used  [1]. 

In  a  separate  study,  under  a  research  agreement  on  flight 
path  reconstruction,  workers  at  the  University  of  Newcastle 
investigated  methods  for  the  rejection  of  deterministic  sinusoidal 
disturbances.  An  optimal  Kalman  filtering  approach  was  found  to 
be  impractical  [2],  However,  additional  research  showed  that  the 
optimal  filter  could  be  approximated,  to  a  high  degree  of 
accuracy,  by  a  digital  filter  which  was  both  simple  and  practical 
to  use  [3,  4,  5]. 

A  simple  program  presented  in  Ref.  3  has  been  substantially 
modified  to  adapt  it  for  use  as  part  of  the  flight  reduction 
procedures  in  the  ABS-RW  group. 

In  Section  2,  the  theoretical  derivation  of  the  filter  is 
summarised  briefly,  while  Section  3  discusses  practical  issues 
involved  with  the  application  the  filter  to  flight  data.  Finally, 
in  Section  4,  the  performance  of  the  new  filter  is  compared  with 
that  of  a  comparable  Butterworth  filter. 
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2.  DERIVATION  OF  THE  DIGITAL  NOTCH  FILTER 

The  aim  was  to  design  a  filter  which  would  remove  a 
sinusoidal  disturbance  from  a  measured  signal,  ym(t),  to  produce  a 
filtered  signal,  yf(t).  In  the  limit,  it  was  desired  that  the 
filtered  signal  should  approach  the  original  signal  y(t)  (see 
Fig.1). 


The  detailed  design  and  the  underlying  theory  of  this  filter 
can  be  found  in  Refs.  2  to  5.  Excerpts  from  the  above  references 
are  reproduced  below. 

It  can  be  shown  that  a  sinusoidal  disturbance,  s,  of 
frequency  u>0>  obeys  the  following  difference  equation  (Ref.  2, 

p2.2). 


D(q'1)s(t)  -  0  (1) 

where  q"1  is  the  delay  operator  and 

D(q-1)  -1-2  cos(ii)oT  )q~,+  q"2  (2) 

Let  the  filter  transfer  function,  F,  be  of  the  form 

F ( q~ 1 )  «  ^4  (3) 

Ds(q  ) 
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1-a.2  cos(w  T)q 
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with 


where  a  Is  a  parameter  such  that  0£  a  SI  and  u 
frequency  being  filtered. 


is  the 


If  the  measured  signal  is  the  sum  of  the  original  signal, 
y(t),  and  the  disturbance,  a(t),  we  have 

yB(t)  -  y(t)*s(t)  (5) 


From  Ref.  1,  page  3*  we  have  the  following  relationship 
between  the  filtered  and  original  signals 

Dg(q-1)  (yf(t)-y(t)}  -  {D(q" 1 )-D3(q“ 1 )}  y(t)  (6) 

Hence,  if  Dg  can  be  chosen  to  be  close  to  D,  while  still 
retaining  stability  in  the  system  as  a  whole,  yf(t)-y(t)  will,  in 
the  limit,  be  small. 

Using  Equations  (M)  and  (6),  this  similarity  between  Dg  and 
D  can  be  achieved  by  selecting  a  value  for  ct  close  to  1  . 

The  final  form  of  the  filter  can  be  expressed  as 


F(q_1) 


G  Blalli 

Vi'  > 


(1-2  cos(uoT)q_1»q"2> 
(1-a2  cos(u)0T)q_1^aV2) 


(7) 


(8) 


where  G  is  a  scaling  factor  to  ensure  that  the  D.C.  gain  of  the 
filter  is  unity  and  is  given  by 


From  Equation  (8),  it  can  be  seen  that  if  the  disturbance 
frequency,  u,  is  equal  to  w0,  exact  cancellation  of  the 
disturbance  will  occur  regardless  of  the  value  of  a  being  used. 
However,  for  u  not  equal  to  <*Qt  the  gain  of  the  filter  is  close  to 
one,  provided  a  is  close  to  one,  so  that  the  signal  passes  with 
minimal  alteration.  Thus  the  filter  resembles  a  notch  at  u 

o 

with  a  a  measure  of  the  notch  width. 


3.  USE  OF  THE  NOTCH  FILTER 

3.1  Discussion  of  Program 

The  original  program  presented  in  Ref. 3  has  been  modified 
and  integrated  into  an  existing  helicopter  data  analysis  program 
[6],  Improvements  were  necessary  to  cope  with  the  cascading  of  a 
large  number  of  filters  and  to  obtain  estimates  of  initial 
conditions.  The  equations  used  in  the  program  are  contained  in 
Section  2.  A  complete  listing  of  the  ARL  program  can  be  found  In 
Appendix  A,  while  Fig. 2  shows  the  flow  of  the  program  during 
execution. 

Execution  of  the  program  requires  the  user  to  type  two 
responses  (the  name  of  the  file  containing  the  data  to  be  filtered 
and  whether  the  filter  characteristics  are  required  in  the  output 
file),  with  the  remainder  of  the  execution  being  controlled  by  an 
input  file  called  'INIT.DAT* .  The  contents  of  'INIT.DAT'  are 
explained  in  the  header  of  the  main  program. 

3.2  Selection  of  Parameter  Values  a,  u 

o 

3.2.1  Background 


Implementation  of  the  digital  notch  filter  is  conceptually 
simple  as  only  two  parameters,  a  and  ,  need  to  be  selected. 
Use  of  the  filter  in  the  theoretical  case,  described  in  Section  2, 
where  u  i3  known  exactly,  simply  requires  the  substitution 
of  a  -  i  and  u  .w  into  the  equations.  This  yields  a  filter 


which  exaotly  cancels  out  the  noise  signal.  However,  In  applying 
the  filter  to  real  data,  a  number  of  practical  considerations  need 
to  be  addressed  as  described  below. 

3.2.2  Choice  of  the  Value  for  a 

The  major  problem  is  that  the  noise  frequencies  in  a  record 
of  data  will  very  rarely,  if  ever,  be  known  exactly,  resulting  in 
an  estimate  for  uQ  being  used  in  the  filter. 

Thus,  if  a  value  of  a  close  to  1  is  used,  the  resulting 
narrow  notch  generated  by  the  filter  may  result  in  the  noise 
signal  being  missed  entirely.  Figures  3,  4,  and  5,  show  the  Bode 
plots  for  the  filter  for  varying  values  of  a,  demonstrating  the 
effect  of  a  on  notch  width.  It  can  be  seen  that  as  a 
approaches  1,  the  notch  width  decreases,  and  distortion  to 
frequencies  other  than  the  cut-out  frequency,  uQ,  decreases. 

In  addition  to  its  effect  on  notch  width,  the  value  of  a 
also  affect3  the  transient  response  characteristics  of  the 
filter.  As  shown  by  sef.  2,  Fig. 2,  and  reproduced  here  as  Fig. 6, 
as  the  value  of  a  tends  to  1  ,  the  roots  of  the  denominator 
polynomial,  Dq ( q  ),  approach  the  boundary  of  the  unit  circle. 
Thi3  results  in  lower  filter  damping.  This  is  illustrated  in 
Fig. 7,  which  3hows  responses  of  the  filter  to  a  step  input  for 
different  values  of  a  .  These  transients  can  be  a  problem  if  the 
initial  conditions  are  offset,  or  if  rapid  changes  occur  in  the 
signal.  In  the  former  case,  the  filter  can  be  started  early,  if 
sufficient  signal  is  available,  and  the  oscillations  allowed  to 
decay  before  the  region  of  interest  is  reached.  However, 
excitation  of  transient  oscillations  during  the  record  are  not  as 
simple  to  remove.  Figure  7(b)  shows  the  response  of  the  filter  to 
the  same  step  input  as  that  used  in  Fig. 7(a),  but  with  the  filter 
run  over  the  data  twice,  in  forward  and  reverse  time  direction,  a 
technnique  used  when  cascading  filters,  as  discussed  below.  The 
figure  shows  that  the  effects  of  ringing  persist  longer  at  the 
higher  values  of  a  .  The  possible  presence  of  such  oscillations 
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should  be  ascertained  when  examining  a  filtered  signal,  and 
allowances  made  if  they  do  exist. 

The  practical  use  of  this  filter  involves  reaching  a 
compromise  in  choosing  the  value  of  a  .  While  sharpening  the 
notch  lessens  the  effect  of  the  filter  on  other  frequencies,  it 
decreases  the  chance  of  filtering  the  signal  if  a  slightly 
erroneous  value  of  is  selected,  ana  also  increases  tne  chance 
of  transient  ringing.  A  method  used  to  partially  offset  the 
effect  of  notch  width  is  to  use  a  series  of  cascaded  filters 
instead  of  a  single  filter.  By  cascading  the  filters,  higher 
values  of  a  are  able  to  be  used,  resulting  in  less  distortion  of 
frequencies  other  than  those  being  filtered.  However,  the 
probability  of  'ringing'  also  increases. 

3.2.3  Cascaded  Filter  Spacing 

Associate  with  tne  cascading  of  a  series  of  filters,  the 
spacing  of  the  filters,  dF,  becomes  an  important  parameter.  A 
large  value  of  dF  (depending  on  a  )  can  allow  intermediate 
frequencies  to  be  passed,  aegrading  tne  quality  of  the  filtered 
signal.  Similarly,  if  tne  spacing  is  too  small,  then  computer 
time  is  being  wasteG  by  repeatedly  filtering  the  same  data. 
Figures  d,  9,  10,  ana  11  show  the  response  of  the  filter  when 

caseaaed  between  two  nominal  frequencies,  3o  and  3-7  Hz. 

3.3  Initial  Conditions 

Errors  in  tne  initial  conditions  lead  to  an  initial 
transient  response  in  the  filtered  data  as  explained  in  the 
previous  section.  In  order  to  estimate  the  initial  conditions,  a 
suoroutine  'INITIAL'  can  be  called.  In  'INITIAL',  tne  filter  is 
initialised  by  setting  the  first  two  filtered  points  equal  to  the 
initial  data  point,  and  the  data  is  cnen  run  through  tne  filter 
aiiowing  any  transients  to  damp  out.  Following  tnis,  the  final 
values  are  used  to  re- ini t iat iae  the  filter  wnien  is  tneri  run  in 
reverse  to  b;ve  tne  required  initial  conditions. 


In  practice,  the  resulting  initial  conditions  depend  on  the 
value  of  (i)Q  chosen  for  the  filter:  thus  it  is  important  to 

choose  uQ  to  match  the  disturbance  frequency  so  as  to  minimise 
any  residual  transients.  If  a  wide  range  of  frequencies  is 
present,  then  the  filter  can  be  re-initialised  at  suitable 
intervals  if  necessary. 


3.*»  Phase  Shift 


As  shown  i3  Figures  3,  J| ,  and  5,  a  phase  shift  can  result 
from  the  U3e  of  a  single  filter.  These  figures  show  that  the 
phase  shift  is  a  function  of  the  value  of  a  used.  As  a  tend3  to 
1 ,  the  region  affected  by  the  shift  in  phase  contracts  to  the 
vicinity  of  the  notch  frequency,  uQ,  and  phase  shifts  can  be 
neglected  at  other  frequencies.  As  outlined  in  Section  3.2.2,  the 
use  of  a  single  filter  is  not  a  practical  method  for  the 
application  under  consideration,  in  contrast  to  cascading.  By 
running  the  cascaded  filters  in  opposite  directions,  any  phase 
shifts  introduced  by  the  filters  tend  to  cancel  each  other. 
Comparisons  of  phase  shift  diagrams  of  Figs.  3,  ^ ,  and  5  with 
those  of  Figs.  8,  9,  10,  and  11  show  the  reduction  in  the  net 
phase  shift  of  the  signal  as  a  result  of  the  cascaded  filtering 
process . 


3.5  Example  using  Spectral  Information 

Because  of  the  ability  of  the  filter  to  remove  specific 
noise  frequencies  from  the  signal  record,  use  can  be  made  of 
spectral  analysis  of  the  signal  to  identify  the  frequencies 
present.  For  example.  Fig. 12  shows  a  vertical  acceleration  signal 
and  Fig. 13  the  noise  power  spectrum  of  the  unfiltered  signal. 
Examination  of  the  noise  power  spectrum  shows  that  there  exists 
considerable  noise  with  frequencies  in  the  region  16  to  19  Hz  and 
in  the  1  to  5  Hz  region.  To  illustrate  the  selective  filtering 
capability,  the  16  to  19  Hz  noise  was  removed  using  filters  with 
a  ,  o.99  and  a  spacing  of  0.05  Hz.  The  resulting  output  signal  is 
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shown  in  Fig. I1!.  The  power  spectrum  (Fig. 15)  of  the  filtered 
signal  shows  the  removal  of  noise  with  frequencies  in  the  1 6  to  19 
Hz  region,  while  other  frequencies  are  unaffected. 

4.  COMPARISON  BETWEEN  NOTCH  AND  BUTTERWORTH  FILTERS 

4.1  Overview  of  the  Butterworth  Filter 

The  purpose  of  developing  the  new  notch  filter  wa3  to 
produce  a  filter  which  performed  better,  when  applied  to  actual 
flight  data,  than  the  Butterworth  filters  currently  used.  The 
characteristics  of  these  Butterworth  filters  are  shown  in  Ref. 5, 
page  8  and  are  reproduced  below  as  Table  1 . 

The  Bode  plots  of  the  two  filters,  defined  in  Table  1,  are 
shown  as  Fig. 16a  (the  high  cut-off  frequency  filter)  and  Fig. 1 6b 
(the  low  cut-off  frequency  filter).  These  two  figures  show  that 
signal  frequencies  appreciably  higher  than  the  cut-off  frequency 
(e.g.  to  8  hz  and  higher  in  Fig.l6a)  can  still  be  passed,  although 
in  attenuated  form.  This  characteristic  means  that  if  all  signal 
frequencies  above  a  certain  frequency  needed  to  be  removed  from 
tne  Uatu  record,  then  tne  cut-off  frequency,  f  ,  has  to  be  set  at 
some  point  well  below  that  frequency.  This  may  result  in  the 
distortion  of  tne  base  signal. 

phase  shifts  resulting  from  the  use  of  Butterwortn  filters 
are  greater  than  those  incurred  by  using  cascaded  notch  filters. 
However,  this  can  be  compensated  for  ~y  eitner  of  two  methods 
(i)  shift  the  data  by  a  factor,  x  (see  Hef.5)  or  (ii)  filter  the 
signax  from  Doth  ends  of  the  data  record  to  achieve  phase 
cancellation . 

4.2  Bode  Plot  Comparisons 


rt'nen  comparing  tne  Butterwortn  an.:  the  new  notch  filter,  it 
should  oe  remembered  that  tne  notch  filter  is  character i sed  Dy  a 
cut-out  frequency  while  the  Butterworth  is  character ised  ;  y  a  cut- 
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off  frequency.  The  notch  filter  can  be  made  to  emulate  the  low- 
pass  Butterworth  filter  under  consideration  by  cascading  filters 
to  remove  all  signal  frequencies  above  a  specific  value. 

The  cascaded  notch  filter  used  for  the  comparison 
has  a-0.99  with  filters  between  4  and  1 6  Hz  at  a  spacing  of  0.05 
Hz.  (Fig. 17  shows  a  Bode  plot  of  this  cascaded  filter).  A3 
indicated  by  Fig. 17,  the  filter  completely  removes  the  signals 
beyond  4  Hz  while  having  little  effect  on  the  amplitude  of  lower 
signal  frequencies.  Additionally,  the  phase  shift  of  the  filter 
is  very  small  and  can  be  ignored.  Figure  18  superimposes  the  Bode 
plots  of  the  two  filters  and  shows  that  the  notch  filter  has  a 
more  square  drop-off  and  smaller  phase  shift. 

4.3  Comparison  using  Actual  Flight  Data 

The  performances  of  the  Butterworth  and  notch  filters  can  be 
compared  by  U3e  of  the  filters  on  actual  flight  data.  Two  data 
records  are  used  for  the  comparison,  the  first  being  the  vertical 
acceleration  signal  and  the  second,  the  pitch  rate  signal,  both 
from  flight  number  17014.  The  acceleration  signal  contained  sharp 
variations,  while  the  pitch  rate  signal  was  slowly  varying.  These 
records  represent  the  two  extremes  in  signal  type  likely  to  be 
encountered  in  the  flight  measurements. 

Firstly,  in  filtering  the  vertical  acceleration  signal,  the 
notch  filter  used  wa3  divided  up  into  two  parts.  The  first  part 
comprised  a  series  of  cascaded  filters  between  6  and  30  Hz  witn 
a-0.975  and  dF«0.1  Hz  (designed  to  remove  all  high  frequency 
noise),  while  the  second  part  was  a  series  of  filters  between  3 
and  6  Hz  with  a  -0.99  and  dF-0.05  Hz.  The  second  bank  of  filters 
needed  a  higher  value  of  a  to  minimise  distortion  to  the  base 
signal  (which  appears  to  have  some  frequency  components  in  the 
region  of  2.5  to  3  Hz),  while  still  filtering  out  the  3.5  Hz  rotor 
frequency  noise.  The  Butterworth  filters  used  were  those  defined 
by  Table  1.  Figure  19  displays  the  comparison  between  the  two 
filter  types. 


The  distortion  of  the  vertical  acceleration  signal  as  a 
result  of  using  the  low  cut-off  frequency  Butterworth  filter  is 
clearly  evident  and  would  indicate  an  inadequate  filter  choice  for 
this  type  of  signal,  where  sharp  variations  are  present.  Use  of 
the  high  cut-off  frequency  Butterworth  filter  appears  to  be  a 
better  alternative,  as  its  output  closely  follows  the  trends  of 
the  original  signal.  Noticeable  in  this  signal  is  the  presence  of 
residual  noise  with  approximately  3.5  Hz  frequency,  indicating 
that  a  considerable  component  of  the  original  noise  signal  is 
still  being  passed. 


TABLE  1 


Specifications  of  Standard  Digital  Butterworth  Lowpass  Filters 


Quantity 

Filter  1 

(low  cut-off  frequency) 

Filter  2 

(high  cut-off  frequency) 

T(  s) 

1/60 

1/60 

N 

5 

5 

/(Hz) 

3-5 

40 

A  (dB) 

500 

3-0 

ft  (Hz) 

M2 

4-0 

Oe  (rad/s) 

7-0331 292892882D 4- 00 

2  ■  55189033061 18D4-OI 

r(s) 

0-459 

0-127 

E»Fo 

1,  1  -74782241 33541 D +06 

1,  4-5580802280576D4-03 

Ei.Fi 

5,  — 8-07698I4534886D4-06 

5,  -1 -6621499551 730D4-04 

Ei,  Fz 

10,  1-49533895274O2D4-07 

10,  2  •  48776282 1 8574 D  4-  04 

£j.  Ft 

10,  —1  -38624627 391 92 D  4-07 

10,  —  1  ■  899 147 5064559 D  4- 04 

Et.Ft 

5,  6  -4345340695383  D  4-  06 

5,  7  •  3691 389324960D  4-03 

Ei,  Ft 

I,  —  1  •  1 962698 1 761 38  D  4- 06 

1,  -1-I598727628397D4-03 

Use  of  the  notch  filter  as  an  alternative  to  the  high  cut¬ 
off  frequency  Butterworth  filter  seems  a  more  favourable 
alternative,  as  it  has  the  same  overall  quality  in  the  filtered 
signal  as  the  Butterworth  filter,  with  the  added  advantage  of 
having  the  3.5  Hz  noise  removed  from  the  data  record.  However, 
some  residual  oscillations,  in  the  region  of  sharp  variations  in 
the  original  signal,  still  remain.  These  are  probably  due  to  the 
filter  transient  characteristics  and  depend  on  the  choice  of  o, 
as  discussed  in  Section  3.2.2.  More  commonly,  the  notch  filter 
would  be  used  on  slowly  varying  signals  of  which  the  pitch  rate 
signal  from  this  flight  appears  to  be  typical.  Figure  19(b)  shows 
the  comparison  between  the  unfiltered  pitch  rate  signal,  the  notch 
filtered  signal,  and  the  two  Butterworth  filtered  signals.  The 
high  cut-off  frequency  Butterworth  filter  would  be  unsuitable  for 
this  record  type  as  there  is  still  considerable  noise  remaining  in 
the  filtered  data  record.  However,  the  low  cut-off  frequency 
Butterworth  filter,  although  it  has  removed  the  noise  from  the 
data,  it  seems  to  have  introduced  some  low  frequency  oscillations 
into  the  filtered  data  which  are  not  apparent  in  the  original 
data. 

Figure  19  (b)  indicates  that  the  notch  filter  would  be  an 
acceptable  filter  choice,  as  the  noise  signal  has  been  removed, 
while  the  base  signal  remains  uncorrupted.  The  filter  used 
had  a-0.90,  dF«0.10  Hz,  with  notches  cascaded  between  3  and  30 
Hz.  Thi3  figure  shows  that  if  the  signal  is  slowly  varying,  then 
a  cascade  of  notch  filters  with  values  of  a  in  the  region  of 
0.90  to  0.95  can  be  used  with  adequate  results.  Accompanying  the 
low  value  of  a,  large  filter  spacings  (dF  =0.10  Hz),  and  hence 
fewer  filters,  can  be  used  to  filter  a  record.  This  results  in 
some  improvement  in  the  CPU  time  used  when  compared  to  cascaded 
filters  with  higher  values  of  a  for  the  same  data  record. 

4.4  Comparison  of  Execution  Times  for  the  Two  Filter  Types 

As  implemented  here,  the  notch  filter  processes  the  data 
many  times  as  it  filters,  while  the  Butterworth  filter  processes 


the  data  only  once.  Hence,  despite  its  simplicity,  the  run  times 
of  the  notch  filter  can  be  significantly  longer  than  those  of  the 
Butterworth  filter.  In  the  application  described  in  Section  4.3, 
Involving  a  total  of  330  cascaded  filters,  the  time  taken  for  the 
high  cut-off  frequency  Butterworth  filter  was  less  by  a  factor  of 
4.25  than  the  notch  filter  for  the  same  range  of  filtering.  (1.3 
CPU  seconds  compared  with  5.5  CPU  seconds). 

The  run  time  of  the  notch  filter  is  proportional  to  the 
number  of  filters  in  the  cascade.  For  small  variations  of  the 
noise  frequency,  for  example  if  the  noise  is  confined  to,  say  17 
to  18  Hz,  then  the  notch  filter  would  take  considerably  less  time 
than  the  corresponding  Butterworth  filter. 

5.  CONCLUDING  REMARKS 

A  simple,  near  optimal  digital  notch  filter  designed  under  a 
research  agreement  with  the  University  of  Newcastle  has  been 
adopted  for  use  at  ARL  for  the  removal  of  unwanted  noise  from 
helicopter  flight  data.  The  ARL  implementation  of  the  filter 
allows  it  to  be  used  conveniently  in  either  of  two  modes  of 
operation.  If  the  noise  power  spectrum  indicates  the  presence  of 
noi3e  at  specific  frequencies,  then  a  single  filter  or  cascade  of 
filters  can  be  chosen  to  remove  particular  frequencies  of 
concern.  Guidance  has  been  provided  for  the  choice  of  notch  width 
and  filter  spacing. 

If  the  noise  spectrum  covers  a  wide  range  of  frequencies,  or 
else  if  a  power  spectrum  is  not  immediately  available,  then  a 
block  filtering  technique  can  be  used  to  remove  all  frequencies 
above  a  chosen  value  (low  pass  filter).  A  relatively  wide  notch 
and  wide  filter  spacing  can  be  used  to  reduce  computing  time.  A 
combination  of  the  above  two  operating  modes  can  be  used  if 
required. 

Comparisons  have  been  made  with  two  Butterworth  filters 
currently  used  to  filter  the  flight  data,  and  an  indication  of  the 
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performance  of  the  notch  filter  obtained.  The  comparisons  showed 
that  the  new  filter  has  a  superior  frequency  response,  both  in 
magnitude  and  phase  shift,  when  compared  to  the  existing 
Butterworth  filters.  However,  execution  times  can  be  several 
times  greater  if  a  large  number  of  filters  are  cascaded 
together.  Also,  the  transient  response  characteristics  of  the 
filter  can  cause  problems  with  rapidly  varying  records  if  the 
notch  width  is  very  narrow.  Further  study  could  be  aimed  at 
circumventing  this  potential  problem  area. 

In  general,  the  new  filter  is  relatively  easy  to  use  and  to 
program  for  a  computer.  It  provides  a  useful  addition  to  the 
suite  of  programs  used  in  the  ABS-RW  group  for  helicopter  flight 
data  analysis.  Its  simplicity  makes  it  straight  forward  to  use  in 
other  applications. 
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APPENDIX  A  :  Listing  of  the  filtering  program 


PROGRAM  FILTER 


This  program  cancels  sinusoidal  disturbances  out  of  an 
arbitary  signal  in  a  given  frequency  range.  The  filter  is  of 
the  form! 


G*(l-2cos  wTq**-l  ♦  q**-2 )✓< 1-ALPHA. 2cos  wTq**- 1 «ALPHA**2 .q**-2 ) 
Where  s 


T 

w  : 

FRE  : 
ALPHA: 


Sampl jng  pgriod. 


Lowest  frequency  component  in  the  sinusoidal  disturbance 


sensitivity  to  inaccuracies  in  the 


CO  -  1)  F.  _ 
value  of  FRE. 

Typically  <  0.95  -  0.99> 

A  coefficient  such  that  the  filter's  gain  is  1.0 


Other  variables  in  the  program  are 


N  : 
NF  : 
FRE  s 
FREXPi 


OF 

NPAS 


Length  of  the  data 

Number  of  filters  in  cascade 

Frequency  used  in  the  design  of  the  first  filter  in  Hz 
The  expected  value  of  the  noise's  frequency  in  Hertz 
Used  in  finding  the  initial  conditions 
Increment  in  frequency  from  one  filter  to  the  next 
Typically  <  0.1Hz  -  l.QHz  > 

Number  of  passes  through  the  data  used  to  find 

the  initial  condition  for  the  filter  output.  This  is 

done  for  each  one  of  the  filters. 


Files  used  in  the  program  i 


" INFILE" 
CLEAN. OAT 
INIT.OAT 


Contains  the  UNfiltered  signal. 

Contains  the  filtered  signal. 

Contains  the  data  for  the  different  variables 
used  in  the  program. 


LAST  UPDATE 


24th  June  1965 


IMPLICIT  DOUBLE  PRECISION  C  A-H,  0-2) 

DATA  PI/'S.  141592654/ 

DOUBLE  PRECISION  Y(3),  UERC600Q),  AUXC6000) 
CHARACTER*500  INFILE 
CHARACTER* 1  ANS 
WRITE  (5.2) 


DIGITAL 


NOTCH 


FORMAT  <$(✓). 

1 

2  '  What  is  the  name  of  the  input  file?1//) 
READ  (5,3)  LENGTH,  INFILE 


FILTERING' 


FORMAT  l 
OPEN  (UNI 
OPEN  (UNIT-30;fILE' 


Q,  A500 

T-20 ,FILE“ INF ILE ( 1 : LENGTH ) .STATUS- 


I'll  I  LL  V  1  •  LLItU  I  II  /  >  J  m  I  UJ~ 

ur  i.ii  sunt  i  -✓  w  if  a  uw  INIT.OAT' , STATUS- ' OLD ‘ ) 
OPEN  ( UN  I T-2 , F l LE« • CLEAN . DAT ' , STATUS- ' UNKNOWN 


OLD* ,ERR- 1 ) 
) 


REAO  (30,*)  FRE,  FREXP,  ALPHA,  T,  OF 
READ  (30,*)  NPAS,  N,  NF 


WRITE  (5,5) 

FORMAT  <.//'  Output  the  filters  characteristics 
READ  (5,6)  ANS 
FORMAT  (All 

IF  (  ANS  .EQ.  'Y'  .OR.  ANS  .EQ.  'y'  )  THEN 

WRITE  (2,7)  INFILE  < 1 i LENGTH) 


?  '  *  ) 


FORMAT  (*  Input  filename  :',A20) 
WRITE  (2,6)  FRE.  FREXP,  ALP^A,  ’ 


FORMAT  ( 


,8)  FRE.  FREXP,  ALPHA,  T.  DF,  NPAS,  N,  NF 
*  '  .  '  FR£  -  '  ,Fl5 . 8  '  F^EXP-'  ,  F 15 . 8  ,  ^  , 


END  IF 


ALPHA- ' ,Fl5 . 8 ,/, ' 'T  6  1  ;F15.8,/, 
DF  -',F15.8,/,'  NPAS  - '  ,  15  ,  /  ,  N 


-'  ,15,- 


NF 


Reed  in  the  data  to  be  filtered  from  tape  20 
DO  10  I  -  1,N 

10  READ  (20,2030)  2,  UER(I) 

Check  to  see  whether  the  data  needs  to  have  the 
correct  initial  conditions  to  be  found. 


IF 


(  NPAS  .GT.  0  )  THEN 

CALL  INITIAL  (FREXP .ALPHA , T , NPAS , N , UER , V ) 


Y  ( 3  ) 

Y  ( 2  ) 
ELSE 

Y( 3  ) 

Y  ( 2  > 
END  IF 


Y  ( 1 ) 
Y(2) 


UER ( 1 ) 
UER( 1 > 


onnn  on  nn  n  nnon  on  no  non  non  onnno  noon  n  noon  non'  onnn  nnonoo  non 


This  taction  reprtttntt  the  main  loop  of  the  program. 

FRE  -  FRE  -  OF 
00  2000  11  -  1,  NF 
FRE  -  FRE  ♦  OF 

WRITE  (5 ,20)FRE 

20  FORMAT  ( 'Filter  frequency  :  '  ,F 10.9) 

Calculate  the  value  of  2.Cos(Wk.T)  and  the  filter's  O.C.  gam  'G* 
A  -  2.0  e  DCOS(2.0*Pl*FRE*T) 

G  -  7  1.0  -  ALPHA*A  ♦  ALPHA*ALPHA  )  s  <  2 . 0  -  A  > 

Perfor#  the  filtering  task  for  tha  remainder  of  the 
date  points 

00  100  I  -  3.  N 

E  -  UER(I)  -  A*UER( I  -  1 )  ♦  UERU-2) 

Y(l)  -  ALPHA*A*Y ( 2  )  -  ALPHA* ALPHA* Y(3 )  ♦  E*G 

Over-write  the  old  value  of  Y(3)  to  the  data 

UERU-2)  -  Y(3 ) 

Shift  tha  data  points  Y<3)  and  Y(2)  to  the  next 
data  point 

Y(3  >  -  Y<2 ) 

Y(2 )  -  Y<1) 

100  CONTINUE 

Write  out  the  two  filtered  end  points  to  the  variable 
array  UER 

UER(N-l)  -  Y(3 ) 

UER(N)  -  Y(2> 


*  V  V  * 

.vs  ■  • 


JNV/  /  V 

aVA.SNj 


v-y-\;.\v 

'n**.'  •.*  s'' 


/  v  ;  ;  : 


Modify  the  values  for  Y(3)  and  Y(2)  so  that  they  contain 
the  correct  values  for  the  reverse  sweep  of  the  data 

Y(  2 )  -  Y(3> 

Y<3 )  -  YU) 

Load  into  the  array  AUX  the  reversed  contents  of  the 
array  UER  to  enable  the  backwards  pass  through  the  data 

DO  200  I  -  1,  N 

200  AUXCI)  -  UER(NU-I) 

Copy  the  contents  of  the  array  AUX  into  the  array  UER 

00  300  1  -  1,  N 
>00  UER(I)  -  AUX  <I) 

2000  CONTINUE 


Check  to  see  wHether  there  exists  a  condition  that 
would  result  in  tha  reversal  of  the  data  and  if  so 
ensure  that  the  data  is  written  out  in  the  correct  form 


/  V  V  V 


IREU  •  2  •  INT  C  FLOAT  (NF)  /  2.0  ) 
IF  <  IREU  .£Q.  NF  )  GOTO  2015 

DO  2005  0J  -  1,N 

2005  AUX  (JJ)  -  UER  <  N  ♦  1  -  JO  ) 
DO  2010  JJ  -  1 ,N 
2010  UER  (  JJ  )  -  AUX  (  JJ  ) 


-.vo 

■•■WV 

\-VvV 


2015  CONTINUE 


DO  2020  !  *  1  ,N 

TIME  -  T  •  FLOAT  (1-1) 

2020  WRITE  (2,2030)  TIME,  UER( I ) 

Call  routine  which  calculates  the  CPU  time  used 
in  tha  running  of  the  program. 

CALL  CPTIME 

*2030  FORMAT (2X,2(F15.?,2X) ) 

END  ’ 


»  V  •  .  *, 

vS:S::;;- 

->>>>', 

vVvy 


•  v- 


non  onn  nnnn  n  non  n  nnn  nnnn  nnnn  nnnnn  nnnnnnnnn 


A  A  A  A  ”■ 


SUBROUTINE  INITIAL  I FRE , ALPHA , T , NPAS , H , UAL . Y ) 

»##**######ee#«e#ee*e*#e#*e#*e#e*e*#*»e* *i*i***f**»» •••••••• ••«»••••« 

The  purpose  of  this  subroutine  is  to  find  the  initial 
conditions  for  the  dste  to  enable  a  single  pass  filtering 
to  be  performed. 

Last  update  :  20th  June  1989 

IMPLICIT  DOUBLE  PRECISION  C  A-H,  0-2  ) 

DIMENSION  Y<3>.  UAL <6 0  00  )  ,  AUX(6000> 

DATA  PI  /3. 1415$2654x 

The  following  is  used  to  ascertain  whether  NPAS  is 
odd  or  even  since  for  correct  operation  of  this  subroutine 
NPAS  MUST  BE  ODD  •*•*♦* 

IREU  -  2  *  1NT  (  FLOAT (NPAS) '2.0  ) 

IF  (  IREU  .EQ.  NPAS  )  THEN 
NPAS  -  NPAS  ♦  1 
END  IF 

Initialise  the  values  for  Y(3)  and  Y(2)  to  the  first  value 
in  the  data  to  be  filtered 

Y<2)  -  UAL ( 1 ) 

Y<3)  -  UAL ( 1 ) 

Calculate  the  value  for  2.CosUJk.T)  and  the 
filter's  D.C.  gain  value  'G' 

A  -  2.0  *  DCOS(2. 0*PI*FRE*T) 

G  -  <  1.0  -  ALPHA*A  ♦  ALPHA*ALPHA  >  /  l  2.0  -  A  > 

This  section  represents  the  main  loop  of  the  subroutine 
10  CONTINUE 

DO  100  I  -  3.  N 

E  -  UAL(l)  -  A*UAL( 1-1 )  »  UAL( I -2 ) 

Y<  1 )  -  ALPHA*A*Y<2)  -  ALPHA* ALPHA*Y ( 3 )  ♦  E*G 

Shift  the  data  points  Y<2>  and  Y<3)  for  the  next  calculation 

YC3)  -  Y<  2 ) 

Y(2 )  -  Y( 1 ) 

100  CONTINUE 

Load  into  array  AUX  the  reversed  contents  of  the  array 
UAL  to  enable  the  backwards  pass  to  be  made 

DO  200  I  -  1,  N 
200  AUX  < I )  -  UAL(N*1- 1 ) 

Copy  the  contents  of  the  array  AUX  into  the  array  UAL 

00  300  I  -  1,  N 
300  UAL(I)  *  AUX(I) 


Decrease  the  value  of  NPAS  by  1 

NPAS  -  NPAS  -  1 
YC2)  -  Y(3 ) 

Y<  3 )  •  Y ( 1 ) 

IF  (  NPAS  .GE.  0  )  GOTO  10 

RETURN 

END 
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Figure  1 


General  Overview  of  the  Filtering  Process 
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Figure  2 


Flow  Diagram  for  the  Notch  Filter  Program 
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Figure  4 


Frequency  Response  of  the  Notch  Filter  for 
ALPHA  =  0.850 
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Figure  8  Configuration  of  the  Roots  of  the  Notch  Filter 
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Figure  7(a)  Comparison  of  the  Step  Response  of  the 

Notch  Filter  for  Various  Values  for  ALPHA 
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Figure  7(b) 


Comparison  of  the  TVansient  Responses  of 

the  Notch  Filter  to  a  Step  Input  for  Various  Values 

of  ALPHA 


'  V’V’  ■.-’jTT^ir 


,  ^ 


>1 


■vy:^ 

“IV* 


m 


'  vV  . 


M.PHA  •  •  .»* 

riLTCRS  (CTUCCN  3.3,  3.7  Hi 
or  •  Hz 


RLPMR  •  *.999 

FILTERS  IETUEEN  3.3,  3.7  Hz 
DF  -  0 .05  Hz 


•r «.  w.  r»  r.  trm  r, 


Noubimoob  ibouaa/s  ibmaiao 


Signal 


v. 


•Vv- '.*>  V 

•V  ‘siA 
^.v>.Vv 


.*  /  -■  ,■ 
V'js 
vAV^-S\ 

La  •  m  * 


% 

V  V  V  V  " 

vJvv 

•'.  »  V  V  V 


I. 


Figure  16(a)  Frequency  Response  of  the  High  Cut-Off 
Frequency  Butterworth  Filter 
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